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Abstract 

In this paper, we have considered a spatially flat FRW universe filled with pressureless 
matter and dark energy. We have considered a phenomenological parametrization of the 
deceleration parameter q(z) and from this we have reconstructed the equation of state for 
dark energy u^(z). This divergence free parametrization of the deceleration parameter is 
inspired from one of the most popular parametrization of the dark energy equation of state 
given by Barboza and Alcaniz [30]. Using the combination of datasets (SN la + Hubble + 
BAO/CMB), we have constrained the transition redshift zt (at which the universe switches 
from a decelerating to an accelerating phase) and have found the best fit value of zt- We have 
also compared the reconstructed results of q(z) and uj^z) and have found that the results 
are compatible with a ACDM universe if we consider SN la + Hubble data but inclusion 
of BAO/CMB data makes q(z) and oj^(z) incompatible with ACDM model. The potential 
term for the present toy model is found to be functionally similar to a Higgs potential. 
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1 Introduction 

The discovery of the late-time cosmic acceleration [1, 2] opened up a new field of research in mod¬ 
ern cosmology. A number of theoretical models have been constructed to explain this accelerated 
phenomenon. Most of them are based either on some modification of the Einstein-Hilbert action 
[3, 4] or the existence of new kind of exotic fields in nature, dubbed as “dark energy” (DE). In this 
paper, we will focus on the second aspect and consider DE as the driving agent for the current 
accelerated expansion of the universe which is considered as a hypothetical energy component 
with a large negative pressure. In the last decade numerous DE models have been explored to 
account for this phenomenon (for review, see references [5, 6, 7, 8, 9]). In spite of those efforts, 
however, the true nature of dark energy still remains a mystery. The most popular and simplest 
cosmological DE model is the ACDM model, which is in good agreement with the recent observa¬ 
tional data. The ACDM model is obtained by introducing a cosmological constant A into general 
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relativity, for which the equation of state parameter cua = — 1. However, it suffers from two 
major problems, namely, fine tuning and cosmological coincidence problems [10, 11]. This moti¬ 
vates theorists to develop some alternative models of DE. Scalar field models, viz, quintessence 
models, are widely used in addressing several issues in cosmology due to their dynamical nature 
[12, 13, 14, 15, 16, 17] for which the equation of state parameter evolves dynamically with time 
contrary to the cosmological constant models. Usually, the potential term for the quintessence 
field dominates over its kinetic part in order to generate sufficient negative pressure and generates 
acceleration. A large variety of quintessence potentials have been considered to serve the purpose. 
But, none of these models have very strong observational evidence. Some excellent reviews on 
this topic can be found in [8, 18] and the references therein. 


From the theoretical viewpoint, in the absence of DE, the universe should be decelerating 
as gravity holds matter together. The existence of an early decelerated expansion phase of the 
universe is also supported by the gravitational instability theory of structure formation and of 
big bang nucleosynthesis. The deceleration parameter, thus, must show a signature flip in order 
to incorporate both the scenarios. In other words, a transition from decelerating phase (q > 0) to 
a late-time accelerating phase (q < 0) is necessary to explain the structure formation as well as 
the current acceleration measurements. Recently, various studies have been performed to analyze 
the kinematics of the universe through phenomenological parametrizations of q(z), in a model 
independent way (see references[19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29]). However, in most of 
these cases, the g-parametrization is valid for z « 1 only, and others diverge at z —> — 1 (see 
references[20, 21, 25, 26, 27, 28]). Recently a new parametrization for DE equation of state c o^(z) 
has been considered which diverges neither in the past nor in future [30, 31], however very few 
such divergenceless parametrizations appear in the literature. The basic characteristics of dy¬ 
namical evolution, both static and dynamic, can be expressed in terms of the Hubble parameter 
H and the deceleration parameter q only. Infact these two parameters enable us to construct 
model-independent kinematics of the cosmological expansion [32], Usually, kinematic approach is 
described by a particular metric theory of gravity. It does not depend on the validity of general 
relativity or any model specific assumptions like the matter-energy content of the universe. Mo¬ 
tivated by this facts, in the present work, we propose a simple two-parameter parametric form 
of q(z) to constrain the evolution behavior of DE such that q(z) does not diverge for the entire 
range z G [— 1 , oo] and we study the expansion history of the universe for this divergence-free 
parametrization. The functional form of g-par arnetrization is similar to the one presented in [30] 
for equation of state parameter cThe properties of this parametrization has been discussed in 
detail in the next section. Based on this consideration, we have then solved the field equations 
for this model and have obtained the expressions for different cosmological parameters, such as 
the equation of state parameter (u^z)). We have also constrained q(z), co^(z) and the potential 
term using the combination of SN la, H(z) and BAO/CMB dataset to investigate the various 
properties of this model. 

This paper is organised as follows. In section 2, we have described the basic theoretical frame¬ 
work for the quintessence model of a spatially flat FRW universe. We have chosen a particular 
form of q(z) to solve the field equations and have reconstructed the equation of state parameter 
oj<t,(z). In section 3, we have considered observational datasets of SN la, H(z) and BAO/CMB and 
have studied the constraints on the various reconstructed model parameters and have summarized 
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the main results of this analysis in section 4. Finally, in the last section, some conclusions are 
presented on the basis of the results obtained in this work. 


2 Field equations and Results 


For a universe having space-time curvature R filled with a scalar field 0 having a potential V (0) 
and normal matter, the action is given by 


S 



V (0) + L m 


( 1 ) 


where L m is the Lagrangian of the background matter which is considered as pressureless perfect 
fluid. Here, all quantities are expressed in units of 8ttG = c = 1. 


For a spatially flat FRW space-time with the metric 

ds 2 = dt 2 — a 2 (t)[dr 2 + r 2 d0 2 + r 2 sin 2 dd(j) 2 ] 


( 2 ) 


the Einstein held equations are obtained as 

3 H 2 = p m + ^0 2 + y(0) (3) 

2H + 3H 2 = -i0 2 + f/(0) (4) 

where an overdot indicates differentiation with respect to time t and p m represents the energy 
density of the background matter component of the universe. We choose a spatially hat geometry, 
which is favoured by the updated results of the cosmic background radiation measurement [33]. 
In a FRW background, the energy density p<f, and pressure density p<p of the scalar held 0 are 


P<fr = P<f> = ~ ^(0)> ^<p = — (5) 

Z Z P(f) 

Here cu^ is the equation of state (EoS) parameter of the scalar held. The conservation equation 
for the scalar held and matter held are 


0 + 3 H (p + 



( 6 ) 


Pm + 3 H Pm — 0 (7) 

Among the four equations (equations (3), (4), (6) and (7)), only three are independent equations 
with four unknown parameters H, p rn , 0 and E(0). So we still have freedom to choose one pa¬ 
rameter to close the above system of equations. 


The deceleration parameter is defined as 


Q = 


a 

TJP 


-(1 + 


H 

JP 


) 


( 8 ) 
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where, q > 0 corresponds to decelerating phase whereas q < 0 indicates accelerating phase of the 
universe. Equation (8) integrates to yield 

H(z) = H 0 ex p ^ 1 i + ^ dz ? ) j (9) 

where, Ho is the Hubble parameter at the present epoch and z = ^ — 1 is the redshift. Clearly, if 
the functional form of q(z) is known, then one can obtain information regarding the evolution of 
the Hubble parameter. 

In principle, we can parametrize the deceleration parameter as, 

Q(z) = Qnfn(z) ( 10 ) 

n =0 

where q n ' s are some real numbers and f n (z) are functions of redshift z. This includes constant 
q model, when fo(z) = 1 and f n (z) = 0 (n > 1). It deserves mention that a large number of 
parametrizations of q(z) have been considered in literature [19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 
29],but the exact functional form of q(z) that will describe the whole evolution history of the 
universe is still being searched for. 

In this present work, in order to close the system of equations, we have considered a new 
parametric form of q(z) as 

q(z) = qo + gi "! 1 , + ^ ( n ) 

1 + z z 

In the above parametrization, qo represents the present value of q(z) and the second term repre¬ 
sents the variation of the deceleration parameter with respect to z. The deceleration parameter 
reduces to q(z) = qo + qi , when z » 1 (at high redshift). One can thus obtain a radiation 
dominated universe by suitably choosing values of qo and q\ for z » 1 limit. However for small z, 
the model represents dark energy behaviour. It should also be noted that the above parametriza¬ 
tion generalizes to the well known linear expansion for q(z) = q 0 + qiZ [22], when z « 1 (at 
low redshift). The advantage of this parametrization is that it provides finite value of q in the 
entire range, z G [—1, oo] and is thus valid for the entire evolution history of the universe. So, 
we can use the parametrization for studying the future evolution of the universe also. It is worth 
noting here that the above parametric form for q(z) was inspired from one of the most popular 
divergence-free parametrization of the dark energy equation of state [30]. We have considered the 
functional form of q(z) containing two terms only in order to make the model simple and elegant. 
Functional forms of q(z) containing suitable higher order terms may provide more information 
regarding the evolution history of the universe but will involve complicated algebra and lacks 
simplicity. But these non-trivial parametrizations always opens up possibilities for future studies 
regarding the nature of DE. We can now easily constrain the two parameters qo and q\ by using 
the available observational data. 

Inserting equation (11) into the equation (9), we obtained the expression for the Hubble 
parameter (in terms of z) as 

H(z) = H 0 (l +z) (1+qo \l +z 2 )^ (12) 

For the sake of completeness, solving (7), we obtain 

PmiyZ) Pmo(l T z') 


(13) 
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Here p m o is an integration constant and represents the present value of the matter field density. 
From equations (3), (12) and (13), we obtained the expression for energy density of the scalar 
held as 


P*(z) = 3tf 2 (l + £) 2(1+9o) (l + -T - 3// 0 2 H m0 (l + z) 3 (14) 


Here h2 m o — fjp is current value of the density parameter of the matter held. 

For this model, the expressions for the density parameters of the matter held (fi m ) and the scalar 
held (fl^,) as well as the equation of state parameter (uj,p(z)) can be easily obtained (in terms of 
z ) as 


£W*) = {0{ = *Wi + ^ (1 " 2,ol (i + Z-)-"' 


(15) 

(16) 


CJ0z) = 


q( z ) 


2 (<?o + _ I) 


3[1 - n m (z)] 3 - 3H m0 (l + Z0~ 2 io)( 1 + z 2 )~n 


(17) 


From equation (17), it can be easily seen: 


1. For z « 1, the EoS parameter c o^z) reduces to 


UJ0Z) 


2g 0 - 1 2qiz 

3 — 3fl m ,o 3 — 3h2 m o 


(18) 


which is similar to the linear redshift parametrization of ui^z) given by u0z) = uq + z, which 
has been used for several cosmological analysis (see references [34, 35]). 


2. For z » 1, we have 


OJ0Z) 


1-2 g 0 - 2 q 1 
3 - 3H m0 ;H- 2( K>- 2 '/i 


(19) 


which is the form of u0z) = (where ca 2 , cj 3 and /3 are arbitrary constants). This 

parametrization is similar to the parametrization of u^z) [36] (for z » 1) for suitably chosen 
values of u> 2 , ca 3 and /3. 


Obviously, the new parametrization of uj^z) or q(z) covers a wide range of other models and 
it also shows a bounded behavior for both high and low redshifts. The simplicity of the func¬ 
tional form of uj<f,(z), however, makes it very attractive and simple to study. 


Now equation (4) can be re-written as 

H = — ~(0 + Pm) 


Substituting H 


— (1 + z)H0t in the above equation, one obtains 


1 

3 Hi 


(1 + z) dH 2 
3Hq d0 


f2 m o(l + z) 3 


( 20 ) 


( 21 ) 
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3 H, 


U 2 = l[qiz(l + ^) 3+290 (l + z 2 ) qi ~ 1 + (1 + g 0 )(l + ^) 2(1+<?o) (l + Z 2 ) qi 


~ ^mo(l + z ) 3 

which immediately gives (by replacing 0 = —(1 + z)H^ 


( 22 ) 


= 00 + 


2 qiz' 


+ 


2(1 + %) 


312 


m0 


(1 + z')(l + z' 2 ) (1 + z') 2 (1 + z') 1+2q °{l +z' 2 ) 


11 


dz' 


(23) 


where 0o is an integration constant. Using equations (3) and (21), then the effective potential for 
the quintessence field can be reconstructed as 




3 m 


Hi 6 Hi dz 


(24) 


Now, substituting H from equation (12) in the above equation, we obtain 


3 H; 


h*V(z) = (1 + z) 2 ( 1+q °\l + z 2 ) q ' - lQ m0 (l + z) s 


(1 + z ) r 


q lZ ( 1 + z) 2(1+w) (l + Z 2 ) qi ~ l + (1 + q 0 ){ 1 + ^)( 1+2 ®)(1 + z 2 ) qi 


(25) 


Therefore, we can reconstruct the kinetic term 0 2 and the potential V(z) (in units of critical 
density p cr i t = 3 Hq) using equations (22) and (25) if the value of Q m0 is given. 

Now, we will apply the same limits used for u^(z) in equation (17) to study the quintessence 
potential in more detail. From equations (23) and (25), it can be easily seen for the following 
limiting cases: 


1. For z « 1, the functional forms of 0 2 and V(z) turn out to be 

1 3 

4>(z) = 0o + -—[2(1 + q 0 ) — 317 m o + 2qiz] 2 

6qi 


and 


V(z) = 3 Hi (- - 


2 llm.o Qo 


- ZH^z 


Combining the above two sets of equations we get, 

U(0) = Uo + rK0-0*)§ 


(26) 


(27) 


(28) 


where U 0 = 77 2 ( 1 —12 m0 ), V — — ^r(3qi)s and 0* = 0o — ^^ +9 ° 3 C? . Note that the above poten¬ 
tial is almost similar to the well-known power-law potential for appropriate choices of Vo, rj and 0*. 


2. Similarly, for z » 1, we have obtained 

1 


3 H 2 


V(z) = ( 


2 _ go _ <h_ 

3 3 3 


z 


2(l-t-<7o+<7l) 




(29) 
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and 


< K z ) = 0o + 


[ 

2(1 + Qo + Qi ) 

3 H m o 

L 

z' 2 

gl (1+250+291) 


2 


dz' 


(30) 


From equations (29) and (30), it is not straightforward to obtain the functional form of V(0). 
However, parametric plot of equations (29) and (30) reveals that in the high limit, the form of 
V (0) is a polynomial in 0 of order four. 

As seen from equation (25), the functional form of V(z) obtained involves a number of model 
parameters and thus depends crucially on the values of these parameters. One can obviously 
choose the model parameters arbitrarily and look at the functional behaviour of the dynamical 
cosmological parameters obtained. These arbitrarily chosen values can then be confronted with 
observational data. But here we do the other way around. We first constraint the various model 
parameters using the available datasets and with the best fit values obtained, we try to reconstruct 
the functional dependence of V(z) or H(0). 


3 Data analysis and method 

We shall fit the theoretical model with the recent observational datasets, namely H(z), SN la, 
BAO and CMB datasets. With the combined datasets, we obtain constraints on the model 
parameters and then try to obtain the functional form of H(0) in the subsequent sections. For 
the sake of completeness, in this section we describe very briefly the techniques followed for 
different datasets. 


3.1 Hubble dataset ( H(z )) 

To constraint cosmological parameter using Hubble dataset, the y 2 function is defined as 


x% = £ l - 

1=1 


h th \ Zi )f 


°H\Zi 


(31) 


where h = ^2 is the normalized Hubble parameter. In above equation subscript “ohs” refers to 
observational quantities and subscript “ th v is for the corresponding theoretical ones. Also, the 
error for normalized Hubble parameter is given by 


a h = h 


/ 



+ 


a H 0 \ 

Ho) 


2 


(32) 


where, on and (?h 0 are the errors in H and H 0 respectively. In this work, we have used the 29 
data points of Hubble parameter measurements [37, 38, 39, 40, 41, 42, 43, 44, 45] in the redshift 
range 0.07 < z < 2.34 (see Table 1 of reference [46]). The current value of the Hubble parameter 
Hq is taken from reference [47]. 


3.2 Type la Supernova dataset (SN la) 

Here we have used the recently released Union2.1 compilation (Suzuki et al. [48]) which contains 
580 data points with redshift ranging from 0.015 to 1.414. The y 2 for SN la is defined as (for 



more details see [49]) 


( 33 ) 


2 _ A B ' 2 

Xsn — A — 

where A, B and C are defined as follows 


(34) 

(35) 


C = T .-2 ( 36 ) 

i =1 u i 

where n° bs represents the observed distance modulus while n th the theoretical one and a, is the 
uncertainty in the distance modulus. 


, W° b "(z>) - /■‘“'MY 

A 2 
i= 1 u i 

n ^ i^n*) -/**(*)] 

B = 1 - -- 


3.3 BAO/CMB dataset 


Next, we have used baryonic acoustic oscillations (BAO) [50, 51, 52] and cosmic microwave 
background (CMB) [53] measurements data to obtain the BAO/CMB constraints on the model 
parameters. The required data points are listed in table 1. For the BAO/CMB dataset, the 
details of methodology for obtaining the constraints on model parameters are given in ref. [54], 
The x 2 function is defined as 


Xbao/cmb 


X T C~ 1 X 


(37) 


with 


X 


/ 


V 


dA{z+) 
Dy (0.106) 

D v < 0.2) 
dA(z*) 
D v ( 0.35) 
dA{z+) 
D v ( 0.44) 
dA(z*) 
D v < 0.6) 
dA (z+) 
D v { 0.73) 


- 30.95 \ 

- 17.55 

- 10.11 
-8.44 
-6.69 
-5.45 ) 


(38) 


where d A (z) = f 0 z is the 
dilation scale and z* ~ 1091 


co-moving angular-diameter distance, Dy(z) = d A (z 


H(z) J 


is the 


is the decoupling time. Also, the inverse covariance matrix C 1 is 


given by 


0.48435 

-0.101383 

-0.164945 

-0.0305703 

-0.097874 

-0.106738 \ 

-0.101383 

3.2882 

-2.45497 

-0.0787898 

-0.252254 

-0.2751 

-0.164945 

-2.45499 

9.55916 

-0.128187 

-0.410404 

-0.447574 

-0.0305703 

-0.0787898 

-0.128187 

2.78728 

-2.75632 

1.16437 

-0.097874 

-0.252254 

-0.410404 

-2.75632 

14.9245 

-7.32441 

-0.106738 

-0.2751 

-0.447574 

1.16437 

-7.32441 

14.5022 / 
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zbao 

0.106 

0.2 

0.35 

0.44 

0.6 

0.73 

dA(z*) 

30.95 ± 1.46 

17.55 ±0.60 

10.11 ±0.37 

8.44 ± 0.67 

6.69 ±0.33 

5.45 ±0.31 

Dv(Zbao) 


Table 1: Values of d^Zbao) f or ^ff eren t va ^ ues °f z bao- 


Dataset 

Qo 

Qi 

y 2 . 

A mm 

Constraints on q 0 , qi 
(within la C.L.) 

H(z) 

-0.82 

0.98 

13.51 

-1.01 <q 0 < -0.64 
0.71 < q x < 1.24 

SN la 

-0.57 

0.70 

562.21 

-0.70 <q 0 < -0.43 
0.27 < qx < 1.14 

SN la ± H(z) 

-0.59 

0.67 

581.18 

-0.67 <q 0 < -0.52 
0.51 <qi< 0.82 

SN la ± H(z) ± BAO/CMB 

-0.5 

0.78 

628.27 

-0.52 <q 0 < -0.47 
0.76 < qx < 0.81 


Table 2: Best fit values of qo and q\ for the present model. 


Finally, the total y 2 for these combined observational datasets is given by 

Xtotal = Xh + Xsn + Xbao/cmb (39) 


4 Results 

In figure 1 we have shown the la (68.3%) and 2 a (95.4%) confidence contours for q$ — qi parameter 
set for different datasets, as mentioned in each panel. The numerical results are displayed in table 
2. In figures 2 and 3, we have shown the evolution of the deceleration parameter and equation of 
state parameter as a function of z respectively. In both the plots, the central thick line is drawn 
for the best fit values of the parameters. The dashed and dotted contours indicate the ltr and 2 a 
confidence levels respectively. 

In figure 2, the green thick line indicates the evolution of q(z) in a standard ACDM model. 
The four panels are for different datasets as indicated in figure 2. It is important to note that 
when one considers SN la + H(z) datasets, the parametrized model is compatible (within 2 a 
limit) with ACDM model, however inclusion of BAO/CMB data makes it incompatible. This 
indicates that BAO/CMB data puts tighter constraints on model parameters. Inclusion of even 
higher order terms or some other form of divergenceless parametrizations in q(z) in equation 
(11) may make better compatibility with ACDM model. However for the present parametrized 
model it is observed that q(z) is not much compatible with ACDM model for SN la + H(z) + 
BAO/CMB data. 
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<?o q o 




<70 


<7o 


Figure 1: Observational constraints on the parameters go one? gi for the parametrization given by 
equation (11). The contours show the la and 2a confidence level from H(z ) (upper left panel), SN 
la (upper right panel), SN la + H(z ) (lower left panel) and SN la + H(z ) + BAO/CMB (lower 
right panel) datasets, for a spatially flat universe. The large dot represents the best fit value of g 0 
and qi. 


It is evident form figure 2 that q undergoes a smooth transition from a decelerating to an ac¬ 
celerating phase of the universe for each dataset. We have also found the observational constraints 
on the transition redshift z t along with their best fit values for different datasets and presented 
them in table 3. The results are found to be consistent with the results obtained independently by 
several authors (see references [21, 22, 23, 25, 29, 55]), which states that the universe at redshift 
z < 1 underwent a phase transition from decelerating to accelerating expansion. However, it is 
interesting to note that the contours obtained with joint analysis of (SN la + H(z) + BAO/CMB) 
dataset put tighter constraints as compared to the constraints obtained from SN la and H(z) data 
sets (see tables 2 & 3). In figure 3, we have plotted 0 J<f,(z) for SN la + H(z) + BAO/CMB datasets 
for different values of f2 m0 . If is evident from the plots that for fl m0 = 0.27 (upper left panel of 
figure 3), the model is not at all compatible with ACDM (c= —1), but as we increase the value 
of Q m o, ACDM model is favoured. However at higher z, the situation is completely different. 
At higher values of z, it is seen that the present parametrized model significantly differs from a 
ACDM model (for which = —1). This dynamical behaviour of cas evident from equation (17) 
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z 


z 


Figure 2: Reconstructed deceleration parameter q for this model using various datasets, as indi¬ 
cated in the each panel. The central thick line (red) represents the best fit curve. The dashed 
and dotted contour represents the la and 2a confidence level respectively. The green thick line 
represents the evolution of q in a flat A CDM with f2 m0 = 0.3 and Gao = 0.7. 


is because of the terms involving higher powers of (1 + z) and is characteristics of the dynamical 
nature of dark energy. This obviously may have its implications on the structure formation of 
the universe. However it is obvious from the plots that at z —> 0, the model approaches ACDM 
limit. It is always nice to have a ACDM limit for a toy model to be relevant with observations, 
but as the dynamical nature of dark energy is still unknown, this deviations from ACDM also 
needs attention. For the combined dataset, the equation of state parameter is constrained to 
be —1.02 < (d$ < —0.96 at the 2cr confidence level (with D m0 = 0.32), whereas the best fit value 
is close to = —0.99 at the present epoch (see lower panel of figure 3). This result is almost 
consistent with the recent observational constraints on u obtained by Wood-Vasey et al. [56] 
and Davis et al. [57] at low z (z < 0.3 or so). 

In ref. [58] the authors have reconstructed the various parameters of scalar held (for example, 
potential term, kinetic term and equation of state parameter) from SN la and BAO datasets using 
non-parametric reconstruction method based on a Gaussian Process representation. They have 
found that the ACDM (c = —1) model is consistent at 2 a level for SN la dataset and at la 
level for BAO dataset. Therefore, such non-parametric techniques may give more robust results 
which are independent of any ad-hoc parametrization of q(z). Here in our case, as discussed 
earlier, because of this particular choice of deceleration parameter, ACDM model is not always 
compatible and may be some different divergence free form of deceleration parameter will satisfy 
the same. 
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Dataset 

Zt 

(best-fit) 

Constraints on z t 
(within la C.L.) 

H(z) 

0.75 

0.71 < zt < 0.81 

SN la 

0.70 

0.58 <z t <l 

SN la + H(z) 

0.82 

0.75 <zt< 0.91 

SN la + H(z) + BAO/CMB 

0.54 

0.51 <zt< 0.57 


Table 3: Best fit values of z t for this model. The constraints on z t are also given. The H(z), SN 
la and BAO/CMB data has been used in carrying out the x 2 analysis. 

SN la + H(z) + BAO/CMB SN la + H(z) + BAO/CMB 




z z 

SN la + H(z) + BAO/CMB 



z 


Figure 3: The reconstructed equation of state for this model using the combined datasets. 
The central thick line is for the best fit values of the parameters g 0 an d <?i arising from the joint 
analysis of (SN la + H[z)+BAO/CMB) dataset. The dashed and dotted contour represents 
the la and 2 a confidence level respectively. In this plot, we have chosen Q m0 = 0.27 (upper left 
panel), Q m0 = 0.3 (upper right panel) and f2 m0 = 0.32 (lower panel). 


In figure 4, we have shown the evolution of the energy densities and the dimensionless density 
parameters for the scalar field (solid line) and matter field (dashed line). The plot is for the best 
fit values of the parameters (g 0 , gi) for the joint analysis of (SN la + H(z)+BAO / CMB) dataset. 

The expression for f(z) obtained in equation (23) is very complicated and it is very difficult to 
find out the functional dependence of f(z). So, we could not express V in terms of 0. However, we 
have used SN la + H(z)+BAO/CMB dataset to reconstruct the form of the effective potential 
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z 
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Figure 4: The left and right panels show the evolution of energy densities (in units of critical 
density) and the behavior of density parameters as a function of redshift z respectively. The 
plots are for the best fit values of the pair (qo, q±) arising from the joint analysis of (SN la + 
H(z)+BAO/CMB) dataset. In both the plots, we have chosen f2 m o = 0.27. 


V(z) as well as the kinetic term (0 2 (z)) of the scalar field 0. The reconstructed behavior of V(z ) 
and <p 2 (z) expressed (in units of critical density) are shown in figure 5. From figure 5, we have 



SN la + H(z) + BAO/CMB 



z 


z 


Figure 5: The reconstructed effective potential V(z) and the kinetic term 0 2 are shown as a 
function of redshift z (in units of critical density). The central thick line is for the best fit values 
of the parameters qo and q\ arising from the joint analysis of (SN la + H(z)+BAO/CMB) 
dataset. The dashed and dotted contour represents the la and 2a confidence level respectively. In 
this plot, we have chosen (2 m0 = 0.27. 

found that the kinetic term 0 2 changes very slowly with z as compared to the potential term 
V at the present epoch. However, we have solved equation (23) numerically and have plotted 
the potential V as a function of 0 in figure 6 for different choices of initial conditions (say 0o) 
and H m o- From figure 6, it is evident that the nature of the behavior of V against 0 remains 
same by small change in the value of H m o- It is also evident from figure 6 that the trajectory 
of V (0) is independent of the initial conditions. It only shifts the plot horizontally. For 0o = 2 
and n m0 = 0.27, a simple numerical evaluation shows that the form of the potential can be well 
approximated as 

H(0) ~ — V70 2 + H 2 0 4 + H 0 (40) 

where Vj = 12.7, V 2 = 1.2 and V 0 = 33.8. These values of V 2 and V 0 has been obtained 
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Figure 6: The effective potential V{4>) (for the potential as given in equation (25)) is plotted 
numerically for the present model with q 0 = —0.5, q\ = 0.78, Q m0 = 0.27 (left panel), 0.28 (middle 
panel) and 0.3 (right panel). The thick, dashed and thin lines are for 0 O = 1,2,3 respectively. 

numerically using the curve fitting technique. It may be noted that we require fine tuning of 
the parameters V\, V 2 and Vo to obtain the exact form of the potential V(<f>). According to the 
standard model of particle physics, the Higgs potential is given by 

r(4) = 1(4> 2 - vr = -^4 2 + + f’ 1 (41) 

where v is the vacuum expectation value of the Higgs held <f> and A represents self-coupling of 
the Higgs potential. Recently, in the framework of a non-minimally coupled theory, this type of 
potential have received a lot of attention in the study of the Higgs inflation models [59, 60, 61]. 
Comparing equations (40) and (41), one can say that the canonical scalar held (f mimics Higgs 
scalar held for the present model for some particular choice of initial condition. 

5 Conclusions 

In cosmology, there are a number of proposals to build an acceptable dark energy model to de¬ 
scribe the early-time and late-time scenarios of our universe. In literature, the dynamical models 
of dark energy offer a better framework to investigate the evolution history of the universe. In 
the framework of a spatially hat FRW universe composed of dark energy and a normal matter 
held, in this present work, we have made an attempt to construct a viable dark energy model 
which shows the desired late-time dynamics of the universe. For this purpose, we have considered 
a simple relation as given in equation (11). Of course, this choice is quite arbitrary and often the 
chosen parametrization leads to possible biases in the determination of properties of a particular 
model parameters. Since we are looking for physically viable model of the universe consistent 
with observations, we make this ansatz in order to close the system of equations. We would like 
to mention here that the g-parametrization adopted in this paper is very simple and valid for the 
entire redshift range. One the other hand, if we consider more than two terms in the expansion 
of the unknown function q(z) (see equation (10)), then it would become very difficult to get firm 
results from data analysis. 

A remarkable feature of this toy model is that the deceleration parameter q undergoes a 
smooth transition from a decelerated to an accelerated phase of expansion in the recent past 
for each dataset (see figure 2). As discussed in the previous section, the values of the transition 











15 


redshift z t are in good agreement with the results obtained by many authors [21, 22, 23, 25, 29, 55]. 

With the g-parametrization, we have reconstructed the equation of state c o^{z) and we have 
shown that the EoS parameter covers a variety of scalar held models (including well-known linear 
parametrization of ay*). In this work, we have found that go < 0 and gi > 0 (within 2 a confidence 
level), which shows the expected behavior of q(z) in accordance with the recent observations. 
Here, (jO^z) does not suffer from the divergency problem like CPL parametrization [62, 63] at 
z = — 1 (see equation (17)) and thus is capable of providing the entire evolution history of the 
universe. 

We have also plotted the density parameters (also energy densities) of the normal matter and 
the scalar held as a function of z. We have shown that the density parameter of the scalar held 
composes the major part of the total density parameter (fi tot = fl m + ) of the universe at 

the recent time and the recent transition from a matter dominated to dark energy dominated era 
occurs at z « 0.53. This result is in good agreement with the observations. 

With the SN la + H(z ) + BAO/CMB data, we have also obtained the evolution of the 
reconstructed potential V(z) for the parametrized deceleration parameter. For the sake of com¬ 
pleteness, we have numerically solved for the potential as a function of 0. The relevant potential 
is found, which mimics a Higgs potential (see equation (40)) for some particular choice of initial 
conditions. This, however, requires more detailed and involved analysis because in that case one 
need not bother about the origin of the scalar held, because the Higgs scalar held has proper 
theoretical background according to the Standard model of Particle Physics and may provide a 
new window to probe the nature of dark energy. 
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